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<N : Abstract 

\ A Coulomb gas representation of dense-matter QCD is derived from a dual 

' transformation of the low-energy effective Lagrangian. The point-like charges 

■ Q = ±1 of the gas are identified with the instantons and anti-instantons 

■ of such topological charges. An instanton repels another instanton with the 
K*" . same force as it attracts an anti-instanton, in contrast to the semiclassical 

o\ 

O 

m . 

O . I. INTRODUCTION AND OVERVIEW 



interaction. 



Instantons have been argued to play an important role in the nonperturbative dynamics 
PhI of QCD. Being classical solutions to the field equations, instantons provide the starting point 
Q_i' for a phenomenologically successful semiclassical approach to QCD Yet the instanton 
^ ■ calculus, as this approach is known, lacks a controUably small parameter which would put 
it on a firmer theoretical ground. The most natural parameter — the diluteness of the 
instanton gas — is numerically small in QCD vacuum, but not to the extent that the 
^ \ interaction between the instantons can be neglected. The instanton "gas" is more similar 
d \ to a liquid than a dilute gas. 

Given its importance for the QCD vacuum, it is interesting to study instanton interaction 
in regimes where the instanton gas can be made arbitrarily dilute. In this paper we perform 
such a study for the regime of very high matter density, where the diluteness is controlled 
by the large ratio of the baryon chemical potential to Aqcd- We find that in this regime 
the instantons (/) and anti-instantons (/) interact with each other via a four- dimensional 
Coulomb potential generated by their topological charges. 

Our approach is different from the standard semiclassical ones. A standard instanton 
calculation starts from the classical instanton gas and proceeds to calculating quantum cor- 
rections. In this approach the // and // interactions start as fundamentally different: at the 
classical level an instanton interacts with anti-instantons, but not with another instanton. In 
contrast, our approach incorporates quantum effects from the very beginning. Our starting 
point is the low-energy effective Lagrangian for high-density QCD with two massless quark 
flavors 0: 

L = f'lidovf - u\d,y^Y] + afi'A' cos(^ - d) . (1) 
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This Lagrangian describes the pseudoscalar isoscalar excitations, (similar to the rj meson), 
and the dependence of ground state energy on the QCD vacuum angle 6. Using a well-known 
formal trick, we rewrite the partition function of the theory (|I]) as that of a Coulomb gas. 
We then notice that 6 plays the role of an external (electrostatic) potential for this Coulomb 
gas. Since, in QCD, the 9 angle is conjugate to the topological charge, we identify the 
constituents of the Coulomb gas with instantons and anti-instantons. In this manner, we 
obtain information about instanton interaction at the quantum level, which is not directly 
available from semiclassical calculations. By comparing the size of the instanton to the 
average inter-particle distance in the gas, we check that the instanton plasma is indeed 
dilute for sufficiently large fi. We also check that the Debye volume contains a large number 
of instantons, which means that the instanton plasma is weakly nonideal. 



II. EFFECTIVE LAGRANGIAN AND THETA DEPENDENCE OF 

HIGH-DENSITY QCD 

Let us briefly recall how Eq. (|l]) is derived (more details can be found in Ref. 0). The 
ground state of high-density QCD with Nf = 2 light quark flavors {u and d) is a color- 
superconducting phase with two condensates of left- and right-handed quarks 



{<lRaAr = ec.,e^'e'^"'Y''- (2) 

Here the structure of the condensates in spin (a,/3), flavor and color (a, 6) indices is 
written explicitly. In the ground state X"^ and have definite length (see below) and are 
parallel. If quarks are massless, then at the perturbative level there is a degeneracy of the 
ground state with respect to the relative U(l) phase between and V^, which means that 
the U(1)a symmetry of the classical QCD Lagrangian is spontaneously broken. The effective 
Lagrangian of the Goldstone boson arising from the breaking of U(1)a has the following form 

L = fiido^r - u\d,^r] . (3) 

Since heavier meson excitations have masses of the order of the superconducting gap A, 
this Lagrangian is applicable only for momenta much smaller than A. For large chemical 
potentials /x ^ Aqcd, the leading perturbative values for / and u are 1^,0]: 



f = -^, u^ = -. (4) 



The field if, canonically normalized, is the field for the t] boson: r] = \/2fLp. Strictly 
speaking, rj is not the lightest mode in this theory. There is an unbroken color SU(2) sector 
which confines at a length scale exponentially larger than the scales we consider in this 
paper. However, the coupling of ip to this sector is suppressed by a power of A// since the 
coupling has the form {j]/ f)FF. Such a sector is absent for the Nf = 3 case briefiy discussed 
in Sec. 0. 

Instantons add a contribution to (^ which is given by [0: 

Vinstif) = -a/i^A^ cos{(f - 6) . (5) 
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This form is uniquely fixed by the anomaly-mediated relation of chiral rotations to 6 ro- 
tations, and by the fact that multi-instanton contributions are suppressed at large fi. The 
calculation of the dimensionless function a(/i) is reviewed below. What is important at this 
moment is that a vanishes in the limit fi ^ oo, which means that the mass of the Goldstone 
boson becomes much less than A at large fi. This is the necessary condition for the effective 
theory (|l|), which contains only one field if, to have a region of validity. 

If the quark masses, rriu and rrid, are non-zero, a new contribution appears in the effective 
Lagrangian (|l]). To lowest order in the masses its dependence on ip is also fixed by symmetry: 

Kiass = -b rriumd cos Lp , (6) 

and the coefficient 6 ~ 1 can be calculated using methods of Ref . . 

The 6 parameter in Eq. (^) is the fundamental parameter of the theory which is equal 
to, or at least very close to, zero in our real world. However, we would like to keep this 
parameter explicit because it plays a crucial role in our further arguments. We note here 
that by minimizing the effective potential Vi^st + Knass with respect to ip one can find the 
theta dependence of the ground state energy. For small rriumd -C a/i^, 

E^^{9) ~ -afi'^A'^ - b niumdA^ cos 9, (7) 

in agreement with the fact that massless fermions suppress the topological susceptibility 
(and the theta dependence). 

The calculation of V^nst(v') can be found in Ref. 0]. The starting point is the instanton- 
induced effective four-fermion interaction [^,|rU|,|I[, 



^inst ^ 

3 

+ 32 



_ 3 

{uR\''uL){dRX''dL) - - {uRa^^\''uL){dRa^^\''dL) 



H.c. . (8) 



By taking the average of Eq. (||) over the superconducting state (0), one finds l^nst, and 
confirms that it is proportional to cos{(p — 9) as in Eq. (|^). Taking the average of (|^) in the 
superconducting ground state, where 

(9) 



we find 



2^271 g 

'4 X2 



^nst(^) = - Jdpno{p)(^-'K^p'^ 12|Xpcos(<^-e), (10) 
where no(p) is the instanton (size) density at finite chemical potential, which is given by 
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Taking the p integration in Eq. (p^, we find 

7 / A„„„ \ 29/3 



sxio^rin^vr^^ 



-) ■ (12) 



Aqcd/ V p 

Thus a — * when fi —>■ oo, so at sufficiently large p the instanton calculations are under 
analytical control. 
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III. COULOMB GAS REPRESENTATION 



The effective low energy dense-QCD Lagrangian (|Tp is the sine-Gordon (SG) Lagrangian. 
Many of the special properties of the SG theory apply. One of these properties is the 
existence of a kink-type solution, corresponding to the domain wall in four dimensions 0. 
Another important property is the admittance of a Coulomb gas (CG) representation for 
the partition function. Such a representation was used previously |jl3|] to argue that, in 
zero-density QCD, the low energy effective Lagrangian with U(1)a anomaly represents the 
dual form of instanton contribution to the partition function. However, QCD at zero density 
is in a nonperturbative regime where no theoretical control is possible |T^. In the present 



paper all calculations are under theoretical control and, therefore, some reliable and precise 
statements can be made. 

The mapping between the SG theory and its CG representation is well known. All we 
need to do is reverse the derivation of SG functional representation of the CG in Ref. [IT^ 
The partition function corresponding to the Lagrangian (|1|) is given byQ 

where we introduced 

A.^. (14) 

u 

Le is the Euclidean space Lagrangian. The imaginary time r is rescaled to bring the kinetic 
term into the Euclidean invariant form in new coordinates with xq = ut. Leaving alone the 
integration over (p{x) for the moment, we expand the last exponent in Eq. (0), represent 
the cosine as a sum of two exponents and perform the binomial expansion: 



>U=0 Q=±l 



M 



The last sum is over all possible sets of M+ positive and M_ negative charges Qa = ±1. 
The last line in Eq. ( [T3| ) is a classical partition function of an ideal gas of M = M+ + M_ 
identical (except for charge) particles of charges -|-1 or —1 placed in an external potential 
given by i{9 — ^{x)). It is easy to see that (for a constant or slowly varying potential) the 
average number of these particle per unit of 4-volume {M)/Vi, i.e., the density, is equal to 
A. Thus making A small one can make the gas arbitrarily dilute. 



^To be precise, the path integral in Eq. (|T3| ) should be understood as an integral over low- 
momentum modes of ip only. The upper limit of the momentum of ip is the ultraviolet cutoff of the 
effective Lagrangian (||), which should be taken as some scale smaller than A. Only tree graphs 
contribute to Z so there is no dependence on the precise value of the cutoff. 
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While 9 can indeed be viewed as an external potential for the gas ([T5|), ^{x) is a dynamical 
variable, since it fluctuates as signified by the path integration in (|l^). For each term in 
(p!5|) the path integral is Gaussian and can be easily taken: 



idT^lloQ- e'^^->^='^^^^'^^"''^ . (16) 



We see that, for a given configuration of charges Qa, —iip{x) is the Coulomb potential 
created by such distribution.^ The function G{x) is the solution of the four- dimensional 
Poisson equation with a point source (the inverse of —d'^): 



G{Xa - Xb) 



Thus we obtain the dual CG representation for the partition function (|1^): 



(17) 



E 



(A/2) 



M 



-J d^x^...j d'xM e-^^E"oQ«e-57^^">*-oQ<>Q^^(— 



(18) 



The two representations of the partition function (p!3|) and (18) are equivalent. 



IV. PHYSICAL INTERPRETATION 

The charges Qa were originally introduced in a rather formal manner so that the QCD 
effective low energy Lagrangian can be written in the dual CG form (|TB|). However, now 
the physical interpretation of these charges becomes clear: since Qnet = Ea Qa is the total 
charge and it appears in the action multiplied be the parameter 6 (see Eq. (0)), one con- 
cludes that Qnct is the total topological charge of a given configuration. Indeed, in QCD the 
6 parameter appears in the Lagrangian only in the combination with the topological charge 
density — i^G^jyG^,y/(327r^). It is also quite obvious that each charge Qa in a given configu- 
ration should be identified with an integer topological charge well localized at the point Xa- 
This, by definition, corresponds to a small instanton positioned at Xa- To corroborate this 
identification we notice that every particle with charge Qa brings along a factor of fugacity 
A ~ a (|n|) which contains the classical one-instanton suppression factor exp(— 871^/(7^ (p)) in 
the density of instantons (pH]). 

The following hierarchy of scales exists in such an instanton ensemble for sufficiently 
large /x. The typical size of the instantons p ~ is much smaller than the short-distance 
cutoff of our effective low-energy theory, A^^. Therefore, in our low-energy description 
they are represented by 6{x — Xa) functions. The average distance between the instantons 



^ One notices that the term a = 6 in the double sum ( [l6| ) is dropped. This is the self-interaction 
of each charge. It would renormalize the fugacity A by a factor exp(— G(0)/(/^m)). This factor 
should be dropped as it represents contribution of very short wavelength fluctuations of Such 
fluctuations have to be cutoff at the scale 1/A. The self-energy of the charges comes from a much 
smaller scale, of order which is already calculated and contained in o. 
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r = A^/^ = (y^/iA) is much larger than both the average size of the instantons and 
the cutoff A^^. The largest scale is the Debye screening length in the Coulomb gas, r/5 = 



(-y/a/2yuA//M) ^ ~ (y/aA) ^. This coincides with the static correlation length of the field, 
which differs from the rest mass of the Goldstone by a factor of u. In short: 



Due to this hierarchy, ensured by large /i/Aqcd, we acquire analytical control. 

It is also quite interesting that, although the starting low-energy effective Lagrangian 
contains only a colorless field ip, we have ended up with a representation of the partition 
function in which objects carrying color (the instantons, their interactions and distributions) 
can be studied. In particular, from the discussions given above, one can immediately deduce 
that // and // interactions are exactly the same up to a sign and are Coulomb-like at large 
distances. 

This looks highly nontrivial since it has long been known that at the semiclassical level 
an instanton interacts only with anti-instantons but not with another instanton carrying a 
topological charge of the same sign. As we demonstrated above it is not true any more at 
the quantum level. Indeed, what we have found is that the interactions between dressed 
(as opposed to than bare) instantons and anti-instantons after one takes into account their 
classical and quantum interactions, after integration over their all possible sizes and color 
orientations, after accounting for the interaction with the background condensates X and 
Y in Eq. (^ (which by itself is nonzero due to the same (anti)instantons as well as due to 
purely perturbative interactions), must become very simple at large distances as explicitly 
described by Eq. (p!8|). One could be really amazed how the problem which looks so com- 
plicated in terms of the original bare (anti)instantons, becomes so simple in terms of the 
dressed (anti)instantons when all integrations over all possible sizes, color orientations and 
interactions with background fields are properly accounted for! 

Such a simplification of the interactions is of course due to the presence of an almost mass- 
less pseudo- Goldstone boson r] which couples to the topological charge. When the instanton 
gas becomes very dilute all semiclassical interactions (due to zero modes) cannot contribute 
much, since they fall off with distance faster then the Coulomb interaction mediated by rj. 
On the other hand, when the instanton density increases (at lower baryon densities), the 
Coulomb interaction becomes more screened and, as the Debye length becomes comparable 
to the inter-instanton distances, we lose analytical control. 

Another interesting observation is that the Gauss law manifests itself in the CG partition 
function by suppressing all configurations with nonzero net charge Qnet- (The weight of such 
configurations is suppressed by exp(— ^Q^^^/L^), where L is the linear size of the system. 
This restricts fluctuations of Qnet to 0{L), negligible compared to normal thermodynamic 
size 0(L^).) In QCD vacuum this corresponds to the suppression of the topological charge 
fluctuations in the chiral limit. Nonzero quark masses give the rj particle a mass, turning 
the long-range Coulomb interaction into a Yukawa one, thus unfreezing fluctuations of Qnet, 
in agreement with a well known result that topological charge susceptibility is proportional 
to quark masses (to the product niumd in our case (0)). 




(size) <^ (cutoff) ^ (distance) ^ (Debye) 
/i-i < A-i < (V^/iA)-i/2 < (v^A)-- 



(19) 
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V. CONCLUSIONS AND OUTLOOK 



In this paper we study a regime of QCD where the instanton dynamics acquires analytical 
control — the high baryon density regime. We showed that the partition function (plSf) of 
the Coulomb gas of instantons is the dual representation of low-energy high-density QCD 
(0) as far as the ground state structure, in particular, its theta dependence, is concerned. 
We discuss the correspondence between these two representations. The most important 
physical result is a rather straightforward identification of the charges Qa from the statistical 
ensemble ( ]I8| ) with the well-known BPST-instantons |T5| , p!6| which have been under intensive 



study since 1970's. We have exphcitly demonstrated that at large distances instantons and 
anti-instantons interact like point-like charged particles. This is in contrast with standard 
semiclassical approximation when II and II interactions are very different from each other. 

A similar correspondence between a formal expression for a statistical ensemble of par- 
ticles and the instanton quarks suspected long ago [0 was conjectured recently [0 based 



on the analysis of the multi-instanton measure. As we have said, no theoretical control was 
possible in Ref. to justify that conjecture. Only in the high-density regime considered 



in the present paper the relation between t] physics and dual representation of the ensemble 
of instantons can be made precise. The properties of the instanton ensemble which we find 
should be of interest for numerical instanton simulations. 

It is possible to generalize our results to the color-flavor-locking (CFL) state of Nf = 3 
QCD ||18[. The U(1)a symmetry is also spontaneously broken in this case. The role of 
the 7] boson is played by the i]' meson, which is also light at high densities [^,^. The 
instanton-induced f]' potential has a form similar to (|^) |]19|,0 : 

^nst iv) = -a' ■ fi^A^ cos(v9 - e) , (20) 

where the evaluation of dimensionless function a' is very much the same as our calculation 
of a in Sec. ^ or Ref. 0]. We need to insert an extra factor m^p into (|TDp and use Nf = 3 
in g{p) in {^A\ ) to find: 

Finally, the same analysis can be also carried in the case of QCD at large isospin density 
fij |2^. The instanton interactions can be also shown to be Coulomb-like and the diluteness 



of the gas is achieved at large /i/. In contrast to finite baryon density case, such a system 
can be simulated on the lattice today. 
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